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ABSTRACT 

We report five new measurements of central black hole masses based on STIS and WFPC2 observations 
with the Hubble Space Telescope and on axisymmetric, three-integral, Schwarzschild orbit-library kinematic 
models. We selected a sample of galaxies within a narrow range in velocity dispersion that cover a range of 
galaxy parameters (including Hubble type and core/power-law surface density profile) where we expected to 
be able to resolve the galaxy's sphere of influence based on the predicted value of the black hole mass from the 
M-o relation. We find masses for the following galaxies: NGC 3585, Mbh = 3.4+o g x 10^ M,,; NGC 3607, 
Mbh = l-2io:4 X 10^ M,v,; NGC 4026, Mbh = 2.1^0:4 x 10^ M,,,; and NGC 5576, Mbh = 1.8+0:4 ^ M,,,, 
all significantly excluding Mbh = 0. For NGC 3945, Mbh = 9+2^ x 10^ M,; , which is significantly below 
predictions from M-o and M-L relations and consistent with Mbh = 0, though the presence of a double bar in 
this galaxy may present problems for our axisymmetric code. 

Subject headings: black hole physics — galaxies: general — galaxies:nuclei — stellar dynamics 



L INTRODUCTION TO BLACK HOLE MASS MEASUREMENTS 

This paper is the latest in a campaign to model the cen- 
tral regions of galaxies to determine masses of putative black 
holes from images and spectra taken with the Hubble Space 
Telescope (HST). The discovery of the presence of a massive 
dark object (probably a black hole) in almost all galaxies hav- 
ing bulges, and the scaling relations found versus host-galaxy 
properties, will stand as one of the important legacies of this 
great observatory. 

Masses were first determined from stellar velocity mea- 
surements made with ground-based telescopes having the best 
possible spatial resolution, together with isotropic kinematic 
models (Dressier & Richstone 1988). When combined with 
the spatial resolution of HST, this method has become the 
standard for black hole mass measurements (e.g., van der 
Marel et al. 1998; Gebhardt et al. 2000b). Black hole masses 
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have also been derived from stellar proper motions in our 
Galaxy (Genzel et al. 2000; Ghez et al. 2005), from mega- 
maser measurements of gas disks around central black holes 
(e.g., Miyoshi et al. 1995), and from gas velocity measure- 
ments (e.g., Barth et al. 2001). Reverberation mapping has 
also been used to find virial products in variable active galac- 
tic nuclei (AGNs; e.g., Peterson et al. 2004). 

Direct dynamical masses are the foundation for all scal- 
ing relations used to infer black hole masses in active galax- 
ies; all measures of black hole mass are derived from the di- 
rect, dynamical measurements. Indirect mass indicators, such 
as AGN line widths, are calibrated to reverberation mapping 
measurements (Bentz et al. 2006), which are themselves nor- 
malized against the direct dynamical measurements (Onken 
et al. 2004). 

A central goal of a companion paper (Giiltekin et al. 2009) 
is the accurate measurement of the intrinsic or cosmic scat- 
ter in the M-o and M-L relationships (Magorrian et al. 1998; 
Gebhardt et al. 2000a). We have thus chosen to augment the 
existing sample of Mbh measurements with new determina- 
tions of Mbh for five galaxies selected to fall within a narrow 
range in velocity dispersion. Along with results from the lit- 
erature, this provides a number of galaxies in a narrow range 
in velocity dispersion large enough that we may probe the in- 
trinsic scatter in the relation without biases incurred by, for 
example, looking only at residuals to power-law fits. 

We present observations of the centers of five early-type 
galaxies in § 2, including HST observations in § 2.3, ground- 
based imaging in § 2.4, and ground-based spectra in § 2.5. We 
report results of dynamical models and black hole masses in 
§ 3, and summarize in § 4. In the Appendix, we provide our 
data tables. 

2. OBSERVATIONS FOR NEW Mbh DETERMINATIONS 

2.1. Observational Sample 

The five galaxies in this study were selected to come from a 
narrow range in velocity dispersion (180 < O < 200 km s^') 
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based on HyperLEDA^^ central velocity dispersion measures 
(Paturel et al. 2003). We chose this range because (1) it in- 
cludes galaxies with Mb ~ —20 mag where both core and 
power-law surface-brightness profiles exist and (2) it includes 
both late- and early-type galaxies. We selected galaxies with 
distances such that the predicted radius of influence was larger 
than 0" 1 . The radius of influence is defined as 



GMbh 
32 (Rinn) '' 



(1) 



where the velocity dispersion a is evaluated at the radius of 
influence. This obviously requires an iterative solution, but it 
converges quickly. The predicted mass comes from the central 
velocity dispersion measurement and the M-a fit of Tremaine 
et al. (2002). 

The sample of new galaxies and their properties are pre- 
sented in Table 1 . Distances are calculated assuming a Hubble 
constant of Hq = 70 km s^ ' Mpc^ ' . We also provide "Nuker 
Law" surface-brightness profile parameters as a function of 
radius given by 



/(r)=2(P-T')/«4 
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which is a broken power-law profile (Lauer et al. 1995). In 
addition to the five galaxies whose black hole masses are re- 
ported in this paper, we observed five others with HST as part 
of the same observing proposal. Two of these (NGC 1374 and 
NGC 7213) had very low signal-to-noise ratio (S/N). The re- 
maining three (NGC 2434, NGC 4382, and NGC 7727) will 
be presented in a future paper. 

2.2. WFPC2 Imaging 

The high-resolution photometry of the central regions of 
the galaxies comes from Wide Field Planetary Camera 2 
(WFPC2) observations on HST using filters F555W (V) and 
F814W (/). The observations, data reduction, and surface- 
brightness profiles (including Nuker profile fits) are detailed 
by Lauer et al. (2005). Surface-brightness profiles are also 
available at the Nuker web page'"*^. 

2.3. STIS Observations and Data Reduction 

Our Space Telescope Imaging Spectrograph (STIS) ob- 
serving strategy and data-reduction methods follow those of 
Pinkney et al. (2003), which may be referred to for details. Ta- 
ble 2 gives the specifications for the STIS observations, which 
used the G750M grating with either a 52" x 0"1 or a 52" x 
0"2 slit along the major axis of each galaxy and the STIS 
1024x 1024 pixel CCD with readout noise of ^ le^ at a gain 
of 1.0. For line-of-sight velocity distribution (LOSVD) fitting 
and for measuring the STIS point-spread function (PSF), we 
used the previously observed stellar spectral templates from 
Pinkney et al. (2003) and Bower et al. (2001), which consist of 
a y = 4.64 mag G8 III star (HR 6770), a V = 5.03 mag K3 III 
star (HR 7576), and a V 3.909 mag KO III star (HR 7615). 
The template stars were scanned across the slit to mimic ex- 
tended sources. 

Most of the STIS setups used an unbinned CCD format 
with a read noise of ~ 1 e^pixel^^. Wavelength range for 
all spectra was 8275-8847 A(Leitherer et al. 2001, pp. 231, 

Available at http://leda.univ-lyonl.fr/. 

See http://www.noao.edu/noao/staff/lauer/wfpc2_profs/. 



234). Reciprocal dispersion measured using our own wave- 
length solutions was 0.554 A pixel^ ' . The distribution of dis- 
persion solutions for a given data set had a a ~ 1.5 x 10^^ 
A pixel^' . The average dispersion given in the Handbook for 
G750M is 0.56 A pixel" ' . We found a comparison line width 
of a = 0.45 A = 17.5 km s"^. Instrumental line widths were 
measured by fitting Gaussians to emission lines on compari- 
son lamp exposures. This gives an estimate of the instrumen- 
tal line width for extended sources. We use approximately 
five lines per exposure, and at least five measurements per 
line. While the comparison lamp exposures were unbinned, 
the galaxy spectra were binned, which increases the measured 
widths by -25% for the 0"1 sUt and by - 3% for the 0"2 slit 
at 8561 A. Leitherer et al. (2001, p. 300) give the follow- 
ing instrumental line widths for point sources: a=13.3, 15.0, 
16.7 km s"^ for the first three G750M setups in Table 2. The 
spatial scale is 0"05597 pixel" ' for G750M at 8561 A. 

The STIS data reduction was done with our own programs 
and FITSIO subroutines (Pence 1998; Pinkney et al. 2003). 
The raw spectra were extracted from the multidimensional 
FITS file, and a constant fit was subtracted from the over- 
scan region to remove the bias level. Because the STIS CCD 
has warm and hot pixels that evolve on timescales of < 1 d, 
the dark-current subtraction needs to be accurate. We used 
the iterative self-dark method described by Pinkney et al. 
(2003). After flat fielding, we vertically shifted the spectra 
to a common dither, combined them, and then rotated. One- 
dimensional spectra were extracted from the final spectro- 
gram using a bi weight combination of rows. A 1 -pixel wide 
binning scheme was used near the galaxy center to optimize 
spatial resolution. We present Gauss-Hermite moments of the 
velocity profiles in the Appendix. 

2.4. Ground-Based Imaging 

CCD images of one of our galaxies, NGC 3945, were ob- 
tained from the MDM 1.3-m McGraw-Hill Telescope. The 
2048 X 2048 pixel CCD named Echelle was used. This chip 
has 0.508" pixel" ' and a readout noise of 2.7 e pixel . 
The conditions were clear but not reliably photometric at all 
times, and the seeing varied between 1"6 and 2"4. The com- 
bined images all had a full width at half-maximum intensity 
(FWHM) of 2"0±0"1. (The surface-brightness profiles at 
small radii are taken from HST data, so the relatively large 
FWHM is not a serious problem.) NGC 3945 was observed in 
/, y , and R filters. Standard CCD reduction tasks within IRAF 
were used to subtract overscan, trim overscan, and divide by 
flat-field frames. The task cosmicrays was used to remove 
cosmic rays because the number of exposures was too small 
for median filtering to work with the V and R bands. The / 
band suffers from interference fringes which did not flatten 
out. The flat fielding was good to 1 % of the sky background 
in R, 2% in /, and 0.6% in V . We present the V-band surface- 
brightness profile for NGC 3945 profile in the Appendix. 

We also obtained wide-field V-band surface-brightness pro- 
files from the following sources in the literature: NGC 3585 
is from Bender et al. (1994); while NGC 3607, NGC 4026, 
and NGC 5576 are from Michard & Marchal (1993). 

2.5. Ground-Based Spectra 

Absorption-line spectra were obtained from three ground- 
based telescopes in order to derive stellar kinematics out to 
large radii. Table 3 describes the spectrograph setups. The 
instrumental resolution, as estimated from widths of compar- 



A Quintet of Black Hole Mass Determinations 



TABLE 1 
Observational Sample 



Galaxy 


Type 


Distance 


Mv 


Profile 


rt 


h 


a 


P 


1 


N3585 


SO 


21.2 ± 1.8 


-22.01 


A 


37.0 


14.72 


1.62 


1.06 


0.31 




E 


19.9 ± 1.6 


-21.56 


n 


70.3 


16.87 


2.06 


1.70 


0.26 


N3945 


SBO 


19.9 ±3.0 


-21.14 


\ 


3.9 


18.62 


0.30 


2.56 


-0.06 


N4026 


SO 


15.6 ±2.0 


-20.32 


\ 


3.0 


15.23 


0.39 


1.78 


0.15 


N5576 


E 


27.1 ± 1.7 


-21.67 


n 


549.2 


17.81 


0.43 


2.73 


0.01 



Note. — Distances are given in Mpc assuming a Hubble constant of Hq = 
70 km s"' Mpc^'. All distances come from surface brightness fluctuations by Tonry et al. 
(2001) except for NGC 3945, which comes from group distance by Faber et al. (1989). The 
uncertainties to distances include random errors only. The third column gives V-band abso- 
lute magnitudes taken from Lauer et al. (2005) and may be converted to V-band luminosities 
via log(Lv/L,r, v) = 0.4(4.83 - My) (see also Verbunt 2008). The fourth column indicates 
surface-brightness profile type: power law (\), core (n), or intermediate (A) as determined by 
Lauer et al. (2005). "Nuker-law" surface-brightness profile parameters are given in columns 5 
through 9 and correspond to Equation 2, where rj; is the break radius in units of pc, 4 is the 
surface brightness at the break radius in units of V magnitudes per square arcsecond, a sets the 
sharpness of the profile break between the outer portion of the profile, which has power-law 
index of (3, and the inner portion of the profile, which has power-law index of y (Lauer et al. 
1995, 2005). 

" NGC 3607 was listed incorrectly as being at a distance of 10.9 Mpc by Lauer et al. (2005). 



TABLE 2 

STIS Long-Slit Spectrograph Configurations 







Slit size 


Exposure 


Name 


Grating 


"x" 


s 


NGC 3585 


G750M 


52x0.1 


12241 


NGC 3607 


G750M 


52x0.2 


26616 


NGC 3945 


G750M 


52x0.2 


22002 


NGC 4026 


G750M 


52x0.1 


9973 


NGC 5576 


G750M 


52x0.1 


7138 



Note. — A summary of the main details of the STIS obser- 
vational set up. Details can be found in the text. 



TABLE 3 
Ground-based Spectrographs 





MDM 


Magellan I 


MageUan II 




2.4-m 


6.5-m 


6.5-m 


Spectrograph 


ModSpec 


B&C 


B&C 


CCD 


Wilbur 


Tek 1 


Marconi I 


Central X 


~ 8500 A 


~ 5175 A 


~5175 A 




Call 


Mgb 


Mgb 


Line widths (a)^ 


i.iA 


0.97A 


0.95 A 




39 kms"' 


56 km s"' 


55 km s"' 


Dispersion (A pixel^') 


0.9 


1.4 


0.78 


Slit length'' 


500" 


70" 


60" 


Seeing (FWHM) 


1.3" 


0.65" 


0.6" 


Slit Width 


0.8" 


0.71" 


0.71" 


Spatial Scale 


0.371"pixer' 


0.44"pixer' 


0.25"pixer' 



" Measured from comparison lamp emission lines.'' As limited by CCD 
format. 

ison lines, was below 60 km s^' in all cases. This allows the 
galaxy line widths (o w 200 km s^') to be easily resolved. 
Using the Modulai" Spectrograph (ModSpec) at MDM Obser- 
vatory with the Wilbur CCD we are able to observe the near- 
infrared Ca II triplet without fringing. At Magellan, however, 
fringing was a problem and so the Mg b spectral range was 
used. Our seeing estimates came from consecutive star ob- 
servations using the same setup. These were confirmed by a 
seeing monitor in the case of Magellan. 

Table 4 summarizes the observations of NGC 3945 and 
NGC 5576. The same, basic observing procedure was used 
at Magellan and MDM. At the beginning of the run, the slit 
width was set to values typical of the seeing, ^ 0"7-l"0. Bias 



frames, continuum lamp flats, twilight sky flats, and compari- 
son lamp spectra were taken before and/or after the night. Cal- 
ibration frames were also taken consecutively with the galax- 
ies. These included comparison lamps, template stars, and 
focus stars. At MDM, the focus frames were created by mov- 
ing the star to new positions along the slit during the pauses 
between 5 and 7 subexposures. These frames allowed us to 
model the spatial distortions in the galaxy exposures more 
precisely. At Magellan, we only had the galaxy peaks and 
single-star exposures to define the "S-distortion." During sub- 
sequent runs, however, we used a flat with a multislit decker 
to see that the S-distortion of a star near the edge of the slit 
would be the same as the S-distortion at the center within one- 
pixel width. We estimate the spatial distortions should be less 
than ~ 0"5 from center to edge. 

For each galaxy, we obtained at least two slit position an- 
gles (P.A.s) to improve spatial coverage. The targets were ob- 
served within |HA| < 2 hr to minimize atmospheric refraction 
and extinction. Multiple galaxy exposures were made at each 
slit position to improve the S/N and to median-filter cosmic 
rays. Some dithering was employed to lessen the impact of 
chip defects. We first obtained exposures at the major axis slit 
position. After two to five exposures, we started the rotated 
exposures. For Magellan I with Tekl, a 600 s exposure gave 
us S/N = 50 per A per 1" wide bin near the Mg b line at the 
center of the galaxy. Only two exposures were required for 
good S/N even at larger radii (where many rows were binned). 
For MDM, a similar (600 s) exposure gave only S/N « 18.4 
per per 1" bin. We therefore used 1200 s exposures (S/N 
w 25.9) at MDM and aimed for more exposures. 

The two-dimensional (2D) spectra were reduced using the 
IRAF tasks found primarily in the ccdred, and twodspec 
packages. Bias subtraction was generally not important, but 
we overscan-corrected and trimmed the CCD frames. Flat- 
fielding was performed using frames constructed out of three 
flats. First, a twilight flat was used to define the illumination 
pattern of sky light along the slit. Second, the small-scale 
structure flat was created by dividing the continuum lamp by 
a smooth fit. Third, the large-scale structure along the dis- 
persion axis was defined in a flat created by fitting ID poly- 
nomials to each row of the continuum lamp flat. These three 
normalized flats were multiplied to obtain the final flat. Since 
the continuum lamp does not have a perfectly flat spectrum. 
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TABLE 4 

Ground-based Observations 



Name 






P.A. 


Exposure 


Type of 


NGC 


Date " 


Telescope 




(s) 


Observation 


3945 


3/16/01 


MDM 2.4 m 





3x1200 


CaT 


3945 


3/16/01 


MDM 2.4 m 


90 


3x1200 


CaT 


3945 


3/17/01 


MDM 2.4 m 


90 


2x1200 


CaT 


3945 


5/10/03 


MDM 2.4 m 





3x1200 


CaT 


3945 


3/15/01 


MDM 1.3 m 




3x300 


/ image 


3945 


3/15/01 


MDM 1.3 m 




2x300 


V image 


3945 


3/15/01 


MDM 1.3 m 




1x300 


R image 


5576 


6/22/01 


Magell 6.5 m 





1x600 


Mgb 


5576 


6/22/01 


Magell 6.5 m 


60 


2x600 


Mgb 


5576 


4/07/03 


Magel2 6.5 m 





2x1200 


Mgb 



" The date of observation given as MM/DD/YY. The observatory and tele- 
scope. The position angle of the slit relative to the major axis of the galaxy. 

Number of exposures x exposure length (s).*^ Type of observation. CaT = 
includes Ca II triplet near 8500 A, Mg b = includes Mg b feature near 5175 

A. 



flat-fielding does not produce an accurate galaxy continuum. 
However, the galaxy spectra are normalized before LOSVDs 
are drawn, so this does not significantly interfere with our 
kinematics. 

After flat fielding, the spectra were wavelength-calibrated 
and corrected for spatial distortions. 

This requires first finding wavelength calibrations as a func- 
tion of position along the slit. We used Ar, Ne, and He com- 
parison lamp exposures to define wavelength as a function of 
position. The fit along dispersion axis was typically a fourth- 
order Legendre polynomial. This provided a root mean square 
(rms) residual of ^0.15 A. As discussed above, the spatial 
axis was rectified using the peaks of stars and galaxies. 

The only remaining steps in reduction were sky subtraction 
and combining exposures. The sky subtraction was performed 
by subtracting a fit to the counts along each cross-dispersion 
row with a low-order polynomial (usually a constant or a line) 
but excluding the central pixels containing significant galaxy 
light. This method worked well for the MDM data with its 
longer slit. Fortunately, for the Magellan I NGC 5576 obser- 
vation, the Moon was down and the sky does not seriously 
affect our line strengths even if it is not subtracted. For the 
Magellan II observations of NGC 5576, other program ob- 
jects with a smaller spatial extent were used to produce a sky 
spectrum. Finally, we averaged exposures using a cosmic- 
ray rejection option. We only combined exposures if they 
were taken on the same night and with the same slit position. 
Any dithering between exposures was removed by shifting all 
galaxy peaks to a common row or column. We present Gauss- 
Hermite moments of the velocity profiles in the Appendix. 

2.6. LOSVDs 

Our modeling of stellar kinematics is done by comparing 
binned LOSVDs of our models to those derived from the 
galaxy spectra. LOSVDs are calculated at the positions in- 
dicated in Figures 1-5. The STIS spectra probe the inner 
and the ground-based spectra probe the outer regions. 
We combine the LOSVDs extracted from both sets of data 
to give us kinematic descriptions of the galaxies from both 
inner and wide-field regions. We deconvolve the observed 
galaxy spectrum using the template spectrum composed from 
the standard stellar spectra. The deconvolution is done with 
the maximum penalized-likelihood method described by Geb- 
hardt et al. (2003) and Pinkney et al. (2003). 
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Fig. 1.— Gauss-Hermite moments of LOSVDs for NGC 3585. Blue 
crosses are Gauss-Hermite moments of LOSVDs from HST STIS data. Also 
plotted are ground-based Gauss-Hermite moments of LOSVDs along the ma- 
jor axis (red diamonds) and minor axis (green triangles) from Fisher (1997). 
Fisher (1997) binned v and a differently from the third and fourth moments. 
We inteipolated and rebinned them consistently. Because the ground and 
/14 moments have been interpolated, their error bars are larger than their scat- 
ter Though Gauss-Hermite moments are not fitted directly in the model- 
ing, the jagged black lines are the resulting Gauss-Hermite fit to the best-fit 
model's LOSVDs from § 3 for the major axis (solid black) and minor axis 
(dashed black). The best-fit model has Mbh = 3.4 x 10^ M,v, and T = 3.5. 
The STIS spectra show a rise in velocity dispersion toward the center and re- 
quire a black hole to match the increased mass-to-light ratio. The red dotted 
line shows the best-fit model for which Mbh = 0, which has T = 4. 1 . 

Ground-based data obtained from the literature in the form 
of Gauss-Hermite moments with associated errors were con- 
verted to LOSVDs using Monte Carlo simulations to estimate 
the uncertainties in the velocity profile bins. The Monte Carlo 
simulations used 10"* realizations for each LOSVD to sam- 
ple the uncertainties. If the Gauss-Hermite moments corre- 
sponded to an unphysical negative value for the LOSVD, we 
assigned it a value of zero with a conservative uncertainty. 
Based on previous experience, we bin the velocity profile into 
13 equal bins that cover the range of velocities seen in the 
given galaxy. Gauss-Hermite data from opposite sides of the 
galaxy were typically averaged (changing the sign of odd mo- 
ments) because our models are axisymmetric. For NGC 5576, 
however, we used the LOSVDs from both sides of the galaxy 
independently, changing the sign of the velocity so as to be 
used with our axisymmetric model. The kinematic data from 
Fisher (1997) binned the higher-order moments (/13 and /14) 
differently from the lower-order moments (V and a). We in- 
terpolated and rebinned them consistently. The smoothness of 
the data suggests that this does not introduce a large system- 
atic error 
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Fig. 2.— Gauss-Hermite moments of LOSVDs for NGC 3607. Symbols 
are as in Figure 1. Ground-based data are from Bender et al. (1994). The 
best-fit model (black lines) has Mbh = 1.25 x 10** M,;, and T = 7.3. The 
best-fit model without a black hole (red, dotted line) has T = 7.5. 

3. DYNAMICAL MODELS TO ESTIMATE Mbh 

In this section we present the results of models to test for the 
presence of a central black hole. We use the three-integral, ax- 
isymmetric Schwarzschild method to make dynamical models 
of the galaxies. The method is explained by Gebhardt et al. 
(2003) and in more detail by Siopis et al. (2008), but we very 
briefly outline it here. 

First, we use the photometric data with the assumption of 
axisymmetry and a given inclination to find the luminosity 
density of the galaxy. We use the luminosity density with a 
given mass-to-light ratio (T) and a given black hole mass to 
calculate the potential. We then calculate the orbits of repre- 
sentative stars in this potential. From this, we determine the 
weights for the set of orbits that best reproduces the surface- 
brightness profile. With the weighted orbit library, we find 
the LOSVD for a given inclination for comparison with the 
observed spectra. 

The results are summarized in Table 5. For all galaxies, 
we modeled several inclination angles, but in all cases none 
was clearly preferred by the models alone, and all found the 
same black hole mass within the stated uncertainties. We 
present the values from our edge-on models below except for 
NGC 3607, whose image indicates that it is face-on. Measure- 
ment errors are la uncertainties. We generally report two val- 
ues for black hole mass and mass-to-light ratio in the text; (1) 
the best-fit-model value which comes from the single model 
with the smallest and (2) the value obtained from marginal- 
izing over the other parameter. These two values are not al- 
ways exactly equal, but the best-fit model is always within la 
of the marginalized value. We put the marginalized values. 



Fig. 3.— Gauss-Hermite moments of LOSVDs for NGC 3945. Symbols 
are as in Figure 1. Ground data are from our MDM observations. Error bars 
for ground-based data are from variations from one side of the galaxy to the 
other, which dominate the total uncertainty for this galaxy. The black jagged 
lines are from the best-fit model, which has Mbh = and T = 7.2. 

which incorporate our uncertainty in the other parameter, in 
Table 5 and use those values for all subsequent calculations. 

We have investigated the reliability of this orbit- 
superposition program in a number of ways. The results from 
the program have been compared to those from the Leiden 
program developed by McDermid et al. (in preparation), giv- 
ing consistent results for NGC 0821 using the different codes 
on the same data and on different data of the same galaxy. 
Siopis et al. (2008) also tested our method by synthesizing a 
distribution function for a galaxy and running them through 
the modeling program. The models were able to recover the 
black hole mass as well as the orbital configuration. 

Gebhardt (2004) demonstrated that a sufficiently large orbit 
library produces consistent results with other orbit libraries of 
similar or larger sizes. He showed that the black hole mass 
is not influenced by the choice of the weight on entropy in 
the solution process, provided it is sufficiently small. Geb- 
hardt et al. (2003) also showed that a set of objects observed 
at HST resolution and ground-based resolution gives consis- 
tent (but with different precisions) black hole masses when 
only the ground-based data are used to construct models. Ko- 
rmendy (2004) showed that estimates of the mass of the black 
hole in M32 using techniques similar to these but different 
in detail give results consistent with the current value (and 
their own error bars) over a 10-fold improvement in spatial 
resolution. Hence, in the context of these models we be- 
lieve the program returns correct estimates of black hole mass 
and mass-to-light ratio, and we adopt these estimates below, 
even when the radius of influence of the resulting black hole, 
^infl = GMBH/o^(^infl), is less than the resolution of the ob- 
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Fig. 4.— Gauss-Hermite moments of LOSVDs for NGC 4026. Symbols 
are as in Figure 1. Ground-based data are from Fisher (1997), for which the 
/13 and /!4 moments have been interpolated. Because of the interpolation, the 
scatter in the data is less than the error bars. The LOSVDs show a sharp 
increase in velocity dispersion toward the center. The jagged lines are from 
the best-fit model, which has Mbh = 2.2 x 10* M; and T = 4.6. The best-fit 
model without a black hole (red dotted line) has T = 5.6. 

servation. 

There are at least three issues that might lead us to report 
black hole masses that are significantly wrong (i.e., outside 
our error bars). First, the models are axisymmetric by con- 
struction, hence significant triaxiality could lead to an error. 
Triaxial models, of course, can reconstruct the correct struc- 
ture (e.g., van den Bosch et al. 2008). Second, the models 
are assumed to have constant stellar mass-to-light ratios ex- 
cept for the central black hole. An admixture of dark mat- 
ter with a spatial distribution different from that of the lumi- 
nous matter could lead us to determine an incorrect black hole 
mass thus requiring dark matter in the model (e.g., Thomas 
et al. 2007). Finally, Houghton et al. (2006) have argued for 
methods of determining the LOSVD that address the limi- 
tation of our method (maximum penalized likelihood) — that 
it produces LOSVDs with correlated errors. Gebhardt et al. 
(2003) did Monte Carlo simulations that give us confidence 
in our ability to estimate the number of independent data in 
our LOSVDs. Nonetheless, our method could be improved. 

3.1. NGC 3585 

NGC 3585 is an edge-on SO galaxy at a distance of 
21.2 Mpc (de Vaucouleurs et al. 1991; Tonry et al. 2001). 
Images of the galaxy show that it is flattened with a nearly 
edge-on dust ring at its center. The LOSVD profile shows a 
roughly constant velocity dispersion from about 1" to ^ 6" 
of O « 200 kms^'. Inside ~ 0"1 the dispersion rises to 
~ 280 km s^^, indicating a Ukely increase of mass-to-light 



Fig. 5.— Gauss-Hermite moments of LOSVDs for NGC 5576. Sym- 
bols are as in Figure 1. The ground-based data are from our Magellan ob- 
servations. The jagged lines are from the best-fit model, which has Mbh = 
1.6x 10*M: andT = 3.6. The best-fit model without a black hole (red dotted 
line) has T = 4.0. 

ratio toward the center For use in analysis of the M-a rela- 
tion, we compute an effective stellar velocity dispersion. 



(3) 



where Rg is the effective radius, /(r) is the surface-brightness 
profile, and V and a are the first and second Gauss-Hermite 
moments of the LOSVD from a slit of width 1". From the 
ground-based velocity profile, we find an effective stellar ve- 
locity dispersion of =213 km s^'. 

The contours in the Mbh-T plane are plotted in Figure 6. 
The best-fitting model is able to reproduce all of the major 
features in the velocity profile, and in order to produce the in- 
crease in velocity dispersion seen toward the center, a black 
hole is required. The velocity profiles for the best-fit models 
are shown in Figure 1. We marginalize over mass-to-light ra- 
tio and take Ax^ = 1 as our lo uncertainty to find a black hole 
mass of Mbh = 3.4log x 10*^ M.,. For Mbh = 0, the marginal- 
ized increases from the minimum by 28.7, which rules out 
the absence of a black hole at very high significance (bet- 
ter than 99.99% confidence). Marginalizing over black hole 
mass, we find a V-band mass-to-light ratio Tv = 3.4 ±0.2. 
The best-fit model with Mbh = 3.4 x 10** M, and T = 3.5 
has = 55.8. The difference between the best-fit model and 
the best-fit model with Mbh = is not obvious in Figure 1, 
which shows the Gauss-Hermite moments of the LOSVDs. 
For this reason, we present a plot of the cumulative difference 
in x^ between the two models in Figure 7. The value for 
is different from the value shown in Table 5, which reports 
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TABLE 5 

Mass Measurements 



Galaxy 






'Wbh.Iow 




high 


Tv 




AX- 


NGC 3585 


213 


3.4 X 10** 


2.8 X 10** 


4.9 X 


10** 


3.4 ±0.2 


55.8 


28.7 


NGC 3607 


229 


1.2 X 10** 


7.9 X 10' 


1.6 X 


10** 


7.5 ±0.3 


69.1 


10.6 


NGC 3945 


192 


9x 10* 


-1.2 X 10' 


2.6 X 


10' 


6.6 ±0.8 


46.2 




NGC 4026 


180 


2.1 X 10** 


1.7 X 10** 


2.8 X 


10** 


4.5 ±0.3 


78.7 


26.2 


NGC 5576 


183 


1.8 X 10** 


1.4 X 10** 


2.1 X 


10** 


3.7 ±0.3 


319.1 


15.5 



Note. — Results from mass modeling. Effective stellar velocity dispersions are given 
in units of km s^', masses are in units of M t , and Ty is units of M: L,^ {, . The black hole 
masses and mass-to-light ratios are the result of marginalizing over the other parameter. 
'Wbh.Iow and Mbh, high are the la confidence limits on the detected black hole mass. The 
final two columns list x^ of the best-fit model, and the difference between the minimum 
in the marginalized %- and at Mbh = 0. 
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Fig. 6. — Mass modeling X' contours for NGC 3585, assuming edge- 
on inclination. Contours are for Ax^ = 1.00, 2.71, 4.00, and 6.63, which 
bracket individual parameter confidence levels of 68.3, 90.0, 95.4, and 99.0%, 
respectively. Contours have been smoothed for plotting. The square shows 
the best-fit model. Dots indicate parameters modeled. The contours were 
smoothed for plotting. The best-fit model has Mbh = 3.4 x 10* Mg and T = 
3.5. Marginalizing over the other parameter we find Mbh = ^■'^-o t ^ 
and T = 3.4 ±0.2. 

the difference in marginalized yj- whereas Figure 7 is the dif- 
ference between two individual models. The best-fit model 
generally differs from the best-fit model without a black hole 
by a larger amount. The cumulative plot shows that most 
of the difference comes from the central '--^ 1". 

3.2. NGC 3607 

NGC 3607 is an elliptical galaxy at a distance of 19.9 Mpc 
(de Vaucouleurs et al. 1991; Tonry et al. 2001). The galaxy 
has a nearly opaque dust disk toward the center (Lauer et al. 
2005), and the image indicates a nearly face-on profile. For 
this galaxy we assumed an inclination angle of 51°, corre- 
sponding to a true axis ratio of 0.4. The dust obscures a large 
part of the bulge, but the nucleus is still visible and is thus 
suitable for modeling. The velocity dispersion is roughly flat 
with radius, but there is an increased rotation in the inner 0"1, 
indicating a dark mass. The effective stellar velocity disper- 
sion Oe = 229 km s^'. The contours in the Mbh-T plane 
are plotted in Figure 8. The velocity profiles for the best-fit 
models are shown in Figure 2. Marginalizing over T, we find 
a black hole mass of Mbh = l-2lo 4 x 10** M,,. For Mbh = 0, 
the marginalized increases from the minimum by 10.6, 
which rules out the absence of a black hole at the 99.9% con- 
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Fig. 7. — Cumulative x^ difference between the best-fit model (Mbh = 
3.4 X 10** M v and T = 3.5) and the best-fit model without a black hole in 
NGC 3585. The difference in x^ is summed from the outermost region to the 
innermost for each axis. Positive Ax^ indicates that the model with a black 
hole is preferred. The total difference in x^ is different from the value listed 
in Table 5 because this difference is between the two individual models rather 
than the marginalized results. Most of the difference comes from the inner 
1". 

fidence level. Marginalizing over black hole mass, we find a 
V-band mass-to-light ratio Tv =7.5±0.3. The best-fit model 
with Mbh = 1 .25 x 10** M ,, and T = 7.3 has = 69. 12. Fig- 
ure 9 shows the cumulative as a function of radius, which 
indicates that most of the difference comes from the central 
~1". 

3.3. NGC 3945 

At a distance of 19.9 Mpc, NGC 3945 is an SBO galaxy 
with a pseudobulge (de Vaucouleurs et al. 1991; Tonry et al. 
2001; Lauer et al. 2007b; Kormendy & Kennicutt 2004). This 
galaxy also contains both an inner bar system and an inner 
disk (e.g., Erwin 2004; Erwin & Sparke 1999). The bars 
may present problems four our axisymmetric code, which we 
discuss in § 4. The velocity profile shows a rise in veloc- 
ity dispersion toward the center from r = 0"2, but this may 
also be interpreted as a slight dip in velocity dispersion at 
r « 0"2-0"3. It has an effective stellar velocity dispersion 
Oe — 192 km s^^ The velocity profiles for the best-fit models 
are shown in Figure 3. 

The contours in the Mbh-T plane are plotted in Fig- 
ure 10. The results of the kinematic modeling show that this 
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Fig. 8. — Mass modeling X" contours for NGC 3607. assuming an in- 
clination of 51°. Contours are the same as in Figure 6. Contours have 
been smoothed for plotting. The best-fit model has Mbh = 1-25 x 10* Mg 
and T = 7.3. Marginalizing over the other parameter, we find Mbh = 
1.2to[l X 10** Mo and r = 7.5 ±0.3. 
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Fig. 9. — Same as Figure 7 but for NGC 3607. Positive values of Ax' 
indicate preference for the best-fit model (Mbh = 125 x 10* Mg and T = 
7.3). Most of the difference comes from the central 1". 

galaxy is consistent with no black hole at its center. Marginal- 
izing over Y, our estimates of the black hole mass for any in- 
clination do not exclude a black hole mass of zero at the la 
level: Mbh = 9^1] x 10'' M , . The 2a upper Hmit is Mbh < 
3.8 X lO'' Mq, and the 3a upper Umit is Mbh < 5. 1 x 10^ M,.. 

Because Mbh = was allowed for this galaxy, we included 
negative black hole masses in our parameter space coverage. 
This let us consider the full extent of the la error distribution 
on the low-mass side. Our model allows a negative black hole 
mass as long as the total mass inside the smallest pericenter of 
the orbit library is positive. In essence, this produces a delta 
function decrement to the mass density. 

The sphere of influence of the black holes of the mass we 
find is below the resolution limit of our data. The black 
hole mass for Rinn/rres — 0.5 is Mbh = D Q^esO^/G = 9.7 x 
10^ M,;;, where D is the distance to the galaxy and Ores = 0"05 
is the spatial resolution limit of the spectra. Such a black hole 
mass, however, is ruled out by our modeling under our as- 
sumptions of constant mass-to-light ratio and axisymmetry, 
as inside 0"1 it would produce an excess velocity disper- 
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Fig. 10. — Mass modeling x^ contours for NGC 3945, assuming edge- 
on inclination. Contours are as in Figure 6. Contours have been smoothed 
for plotting. The best-fit model has Mbh = -2.5 x 10* M,v, and T = 6.8. 
Marginalizing over the other parameter we find Mbh = 9^21 ^ 1^ ' ^■I'l 
r = 6.6 ±0.8. 
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Fig. 11. — Same as Figure 7 but for NGC 3945. Positive values of Ax^ 
indicate preference for the best-fit model (Mbh = ™d T = 8.4) compared 
to the best-fit model with Mbh = 10 x 10** M,-, the mass for a marginally 
resolved sphere of influence. Most of the difference comes from the central 
0"1. 

sion above that observed, especially in the central STIS pixel. 
Marginalizing over black hole mass, we find Yy = 6.6 ± 0.8. 
The model with Mbh = -2.5 x 10*^ M. and Y = 6.8 has 
=45.4. Fi gure 1 1 shows the cumulative as a function 
of radius, which indicates that most of the difference comes 
from the central ^ 0" 1 . 

3.4. NGC 4026 

NGC 4026 is an SO galaxy at a distance of 15.6 Mpc (de 
Vaucouleurs et al. 1991; Tonry et al. 2001). Images show a 
very flattened profile, indicating that edge-on models are ap- 
propriate for this galaxy. There is also a weak, cold stellar disk 
in the center. The spectra show a flat velocity dispersion from 
r « 3"to r « 0"3 of a w 160 km s"^ Inside 0"3, the velocity 
dispersion increases quickly to a = 258 km s^' at the cen- 
ter, a strong indication of increased mass-to-light ratio. The 
effective stellar velocity dispersion = 180 km s^^ 

The contours in the Mbh-Y plane are plotted in Fig- 
ure 12. The velocity profiles for the best-fit models are shown 
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Fig. 12. — Mass modeling x~ contours for NGC 4026, assuming edge- 
on inclination. Contours are as in Figure 6. Contours have been smoothed 
for plotting. The best-fit model has Mbh = 2.2 x 10** M^. and T = 4.6. 
Marginalizing over the other parameter we find Mbh = 2.11q 4 x 10** M,;, 
and r = 4.5 ±0.3. 



Fig. 14. — Mass modeling contours for NGC 5576, assuming edge- 
on inclination. Contours are as in Figure 6. Contours have been smoothed 
for plotting. The best-fit model has Mbh = 1.6 x 10* M,:, and T = 3.6. 
Marginalizing over the other parameter, we find Mbh = l-8^o 4 x 10* Mq 
and r = 3.7 ±0.3. 
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Fig. 13. — Same as Figure 7 but for NGC 4026. Positive values of Ax^ 
indicate preference for the best-fit model (Mbh = 2.2 x 10* M ., and T = 4.6). 
Most of the difference comes from the central 1". 

in Figure 4. Marginalizing over T, we find a black hole mass 
of Mbh 2.lt|S;^ x 10^ Mq. For Mbh = 0, the marginal- 
ized increases from the minimum by 26.2. Such an in- 
crease in Ax^ rules out the absence of a black hole at a con- 
fidence level greater than 99.99%. Marginalizing over black 
hole mass, we find Tv = 4.5 ± 0.3. The best-fit model with 
Mbh = 2.2 x 10^ M,, and T = 4.6 has = 78.7. Figure 13 
shows the cumulative as a function of radius, which indi- 
cates that most of the difference comes from the central ~ 1". 



3.5. NGC 5576 

NGC 5576 is an E3 radio galaxy at a distance of 27.1 Mpc 
(de Vaucouleurs et al. 1991; Tonry et al. 2001). The nucleus 
is offset by ^ 0"04 from the center of the outer isophotes 
(Lauer et al. 2005). Our ground-based spectroscopy quali- 
tatively confirms this, showing that the central region is kine- 
matically separate from the outer regions. The ground-based 
spectroscopy reveals an effective stellar velocity dispersion of 
<3e = 183 km s^'. Unlike the other modeling, for this galaxy 
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Fig. 15. — Same as Figure 7 but for NGC 5576. Positive values of Ax^ 
indicate preference for the best-fit model (Mbh = 1.6x10* M,:, and T = 3.6). 
Most of the difference comes from the central 0"2. 



we did not average both sides of galaxy for the ground-based 
data. Instead, we included data from both sides of the galaxy 
with the sign of the velocity appropriately changed on one 
side. Figure 14 shows the % contours from dynamical mod- 
els. The velocity profiles for the best-fit models are shown in 
Figure 5. Marginalizing over T, we find a black hole mass of 
Mbh = 1-8^04 x 10** M,,. At Mbh = 0, the marginalized 
increases 15.5 above the minimum, indicating that Mbh = 
is ruled out at the 99.99% confidence level. Marginalizing 
over black hole mass, we nominally find Yy = 3.7 ± 0.3. The 
best-fit model with Mbh 1.6 x IO^Mq and Y = 3.6 has 
= 319.3. The value for %^ is much larger than it is for 
the other galaxies because (1) the ground-based data come 
from 13 velocity bins instead of four Gauss-Hermite mo- 
ments, resulting in more constraints, and (2) there are two 
sets of LOSVDs for each ground-based axis measurement: 
one from each side of the galaxy. Figure 15 shows the cu- 
mulative as a function of radius, which indicates that most 
of the difference comes from the central ^ 0"2. 
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4. DISCUSSION AND SUMMARY 

It is worth emphasizing that our modeling gives results for 
massive dark objects at the centers of galaxies, which we call 
"black holes." Though the cumulative evidence in favor of 
these central dark masses as black holes is strong, alternatives 
can only be ruled out for the most highly resolved sources: the 
Galaxy, M31, and NGC 4258 (e.g.. Miller 2006). Thus, while 
a large black hole is the most astrophysically likely explana- 
tion of these central dark objects, they do not strictly have to 
be black holes. 

Our models assume a constant mass-to-light ratio for the 
stellar component of the galaxy. One possible source of sys- 
tematic error may come from the role that dark matter halos 
play. Dark matter halos likely increase the total mass-to-light 
ratio in the outer parts of the galaxy but have less of an impact 
toward the center Thus, by neglecting the dark matter halo, 
we may be overestimating the stellar mass-to-light ratio at the 
center of the galaxy and, consequently, underestimating the 
black hole mass. Future models will incorporate dark matter 
halos. 

Another assumption of ours that may be violated is that of 
axisymmetry. Triaxiality has been addressed in other models 
(van den Bosch et al. 2008). In addition to triaxiality are bars, 
which affect NGC 3945, a double-barred system. While the 
primary bar is not prominent inside of 15", where all of our 
kinematic data come from, box orbits from bars can travel to 
the center. Our modeling code is axisymmetric and simply 
cannot model bars. It is possible that our modeling results, in- 
cluding the mass of the black hole, are skewed by the bars. If 
the bar is aligned mostly along the line of sight, then the line- 
of-sight velocity would be higher than without the bar The 
higher velocities, which would be observed, could be misin- 
terpreted as extra dark mass since the bar orbits would not be 
accounted for in the model. On the other hand, if the bar is 
aligned mostly perpendicularly to the line of sight, the bar or- 
bits would contribute little to the line-of-sight velocity at the 
center, leading to an underestimate in central dark mass. One 
of the bars in NGC 3945 appears to lie in the plane of the sky, 
and one of the bars is at least partially along the line of sight, 
assuming that they are coplanar with the outer disk. Thus it 
is entirely possible that the bars are leading to an incorrect in- 
ference of the black hole mass, but it is not obvious whether 
it is skewed to a high or low value. 

4.1. Anisotropy 

In Figures 16-20, we show the velocity dispersion tensor 
for the best-fit model for each galaxy by plotting the ratio of 
the radial velocity dispersion (o^) to the tangential velocity 
dispersion, defined as = 0.5(Og + o^) so that a,/af = 1 
for an isotropic distribution. Here, is the second moment of 
the azimuthal velocity relative to the systemic velocity rather 
than relative to the mean rotational speed. Uncertainties may 
be estimated from the smoothness of the profiles (Gebhardt 
et al. 2003) to be 0.1 to 0.3. All galaxies are dominated by 
tangential motion at the center. 

NGC 3585 (Figure 16) has an intermediate surface- 
brightness profile and actually shows Or/Or mildly increas- 
ing toward the center along the major axis, but with a steep 
decrease inside of 0"1 along the minor axis. For almost the 
entire range out to 23", 0,-/Or < 1- 

The two core-profile galaxies, NGC 3607 (Figure 17) and 
NGC 5576 (Figure 20), are both ellipticals and show near- 
isotropic distributions in the outer regions, but both galax- 
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Fig. 16. — Shape of the velocity dispersion tensor for NGC 3585 from the 
best-fit model orbit solution. The black line is along the major axis, and the 
red line is along the minor axis. The values for the central part of the galaxy 
are plotted at a radius of 0!'01. The dotted line shows the radial extent of the 
ground-based spectroscopic data. 
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Fig. 17.— Same as Figure 16 but for NGC 3607. 

ies are dominated by tangential motion at the center. For 
NGC 3607, the radial motion drops to almost zero. NGC 3607 
has a relatively strong velocity gradient across the central 0"3, 
and it has a strong drop in the velocity dispersion in the cen- 
ter The only way to reproduce these observables is to have 
complete tangential anisotropy, i.e., no radial orbits. 

The other two galaxies, NGC 3945 (Figure 18) and 
NGC 4026 (Figure 19) are SO galaxies with power-law pro- 
files. For NGC 3945, outside of 10", where the kinematic 
data end, radial anisotropy dominates, but the uncertainties 
are large here (Gebhardt et al. 2003). Inside of 10", the dis- 
persion along the major axis steadily decreases from a roughly 
isotropic value to Or/o, w 0.5. The dispersion along the minor 
axis jumps from tangential to radial at r w 3" and then steadily 
decreases to Or/Ot « 0.7. NGC 4026 shows Or/Ot < 1 almost 
everywhere. 

4.2. Demographics 

We plot the masses found in § 3 against and Ly in Fig- 
ure 2 1 , along with the M-a and M-L relations from Tremaine 
et al. (2002) and Lauer et al. (2007b), respectively. With the 
exception of NGC 3585 in Figure 21(b), all of the black hole 
masses differ from the values predicted by the scaling rela- 
tions by at least la. This shows that our black hole mass 
measurements are precise enough to probe the intrinsic scat- 
ter in these relations. The measurement of the intrinsic scatter 
in these relations is addressed by Giiltekin et al. (2009). If un- 
accounted systematic errors are large, however, the residual 
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Fig. 1 8.— Same as Figure 16 but for NGC 3945. 
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Fig. 19.— Same as Figure 16 but for NGC 4026. 
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Fig. 20.— Same as Figure 16 but for NGC 5576. 

scatter could be due to these. Random errors in distance are 
unlikely to be a large part of this as they are typically 10%, 
which is substantially smaller than the ~ 2 deviation from the 
M-o ridgeline. For these particular galaxies, inclination does 
not appear to make a significant difference in black hole mass. 
Systematic errors from triaxiality or bars, however, are diffi- 
cult to estimate and may contribute. 

If one assumes an intrinsic scatter of 0.3 dex in the M-o re- 
lation (the maximum intrinsic scatter found by Tremaine et al. 
2002) and 0.5 dex in the M-L relation (Lauer et al. 2007a), 
all black hole masses are consistent with the scaling relations 
except for NGC 3945, which is significantly below both rela- 
tions. The black hole masses expected for NGC 3945 from the 



M-o and M-L relations are 1.1 x 10^ M,, and 1.6 x 10^ M j, 
respectively. The 3a upper limit for NGC 3945 is Mbh < 
5.1 X 10^ M j. Hence, while it not possible to rule out the ex- 
istence of a small black hole in NGC 3945, it does not fall on 
the M-o or M-L relations. 

It is interesting to note that the single galaxy in our sample 
that is consistent with having no black hole, NGC 3945, is 
a pseudobulge. Comparing fits to the M-a relation of pseu- 
dobulges with normal bulges and ellipticals, Hu (2008) con- 
cluded that pseudobulges have systematically smaller black 
holes. The small or absent black hole in NGC 3945 is consis- 
tent with those findings. 

4.3. Summary 

We conclude by summarizing the main results of this pa- 
per. We observed five early-type galaxies with high spatial 
resolution kinematics from STIS, which we combined with 
WFPC2 photometry and ground-based observations of pho- 
tometry and kinematics. We modeled these data with three- 
integral, axisymmetric orbit models and found a black hole 
mass consistent with zero and significantly below the M-o 
and M-L relations in one, 

NGC 3945 : Mbh = 9^21 x 10*' Mq, 

though the presence of a double bar in this galaxy may present 
problems for our axisymmetric code. We find evidence for 
central black holes in the remaining four: 

NGC 3585 : Mbh = 3.4+^;^ x 10^ M,„ 
NGC 3607 : Mbh = 1 .2^0:4 x 10^ M,„ 
NGC 4026 : Mbh = 2.1+[J;^ x 10^ Mp, 

and 

NGC 5576 : Mbh = 1 •8^0:4 x 10** M,,,. 

In all of these last four galaxies, the absence of a central dark 
object is ruled out to very high significance. 

We thank Louis Strolger for obtaining CCD images for us 
during his run on the MDM 1.3 m as well as the anonymous 
referee. K.G. thanks Marta Volonteri and Monica Valluri for 
helpful discussions. This work made use of the NASA's As- 
trophysics Data System (ADS), and the NASA/IPAC Extra- 
galactic Database (NED), which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Adminis- 
tration. Financial support was provided by NASA/HST grants 
GO-5999, GO-6587, GO-6633, GO-7468, and GO-9107 from 
the Space Telescope Science Institute, which is operated by 
AURA, Inc., under NASA contract NAS 5-26555. 
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0-„(kmS') Ly{Lgy) 

Fig. 21 . — Plots of black hole masses as functions of the host galaxies' effective stellar velocity dispersion (left) and luminosity (right). The lines are the M-tJ 
relation from Tremaine et al. (2002) (left) and the M-L relation from Lauer et al. (2007b) (right). With the exception of NGC 3945, all masses are found near the 
lines, but they are all inconsistent with the best-fit M-c relation at the la level, and all but NGC 3585 are inconsistent with the M-L line. Given the estimated 
scatter in the relations, most of the black hole masses follow the relation. NGC 3945, however, is significantly below the relations. Even the 3o upper limit to its 
mass (Mbh <5.1 X 10' M,,) is more than a factor of 2 below the expected 1.1 x 10** M ., from M-a and 1 .6 x 10* M ; from M-L. 
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APPENDIX 



This Appendix gives tables of the original data used in this paper. Details are available in § 2. 



A Quintet of Black Hole Mass Determinations 



TABLE 6 
Velocity Prohle for NGC 3585 



Radius 


V 






a 






hi 


h4 


0.00 


-8.1±31, 


.7 


281, 


.2±29, 


.6 


0. 


,030±{).069 


-0.034±0.061 


0.05 


14.3±35, 


.1 


271, 


.0±36, 


.6 


-0. 


125±0.072 


-0.008±0.066 


0.10 


73.2±29, 


.2 


260, 


.9±33, 


.5 


-0. 


156±0.067 


0.044±0.080 


0.18 


132.6±29, 


.7 


193, 


.5±30, 


.6 


-0. 


,057±0.047 


-0.058±0.049 


0.30 


147.7±29, 


.5 


211, 


.2±24, 


.6 


-0. 


181±0.066 


-0.010±0.054 


0.58 


35.0±25, 


.6 


213, 


.3±15, 


.8 


-0. 


,025±0.039 


-0.096±0.021 


1.12 


150.7±25, 


.7 


185, 


.8±22, 


.5 


-0. 


109±0.060 


-0.015±0.064 


-0.05 


-42.1±26, 


.9 


244 


.4±35 


.2 


0. 


106±0.077 


0.094±0.070 


-0.10 


-112.4±29, 


.6 


185 


.2±71 


.8 


0. 


,016±0.115 


-0.020±0.151 


-0.18 


-88.8±28, 


.9 


197, 


.2±21, 


.1 


0. 


,024±0.035 


-0.060±0.036 


-0.30 


-13.3±46, 


.5 


248, 


.6±34 


.4 


0. 


,185±0.071 


-0.045±0.063 


-0.58 


12.3±37, 


.9 


198 


.9±22, 


.2 


0, 


,147±0.054 


-0.038±0.037 


-1.12 


-104.4±32, 


.0 


239, 


.6±27, 


.5 


0. 


,320±0.091 


0.218±0.107 



Note. — Gauss-Hermite moments for velocity profiles derived from STIS 
data. Radii are given in arcsec, first and second moments are given in units of 

km s"'. 



TABLE 7 
Velocity Profile for NGC 3607 



Radius 




V 






G 




hi 


h4 


0.00 


13, 


.2±22, 


.8 


180 


.5±22, 


.1 


0.077±0.068 


-0.007 ±0.063 


-0.05 


-36, 


.1±29, 


,0 


227, 


.3±35, 


.1 


0.119±0.086 


0.069±0.111 


-0.10 


8, 


.9±26, 


.3 


151, 


.4±27, 


.0 


0.000±0.034 


-0.065±0.066 


-0.18 


-59, 


.9±49, 


.4 


216 


.9±39, 


.8 


-0.167±0.096 


-0.021±0.102 


-0.30 


39, 


.1±33, 


.9 


206 


,9±29 


,5 


-0.071±0.091 


-0.040±0.098 


-0.58 


22. 


.8±31, 


.4 


197 


.2±32, 


.3 


-0.003±0.079 


-0.050±0.072 


-1.12 


49 


.2±27, 


,5 


160 


.2±25, 


.7 


-0.011±0.076 


-0.031 ±0.084 


0.05 


-23, 


.3±29, 


.8 


175, 


.1±28, 


.9 


0.074±0.071 


-0.032±0.055 


0.10 


-67, 


.4±34, 


.2 


197, 


.6±44, 


.4 


0.117±0.105 


0.031 ±0.1 20 


0.18 


-71, 


.4±28, 


.1 


169, 


.9±26, 


.0 


0.031 ±0.050 


-0.073±0.036 


0.30 


-21, 


.7±87, 


.7 


251, 


.7±41, 


.6 


0.222±0.178 


-0.002±0.397 


0.58 


-41, 


.8±44, 


.4 


261, 


.8±38, 


.5 


-0.119±0.121 


-0.097±0.121 


1.12 


-27, 


.7±38, 


.5 


197, 


.3±36, 


.3 


-0.005±0.083 


-0.073±0.073 



Note. — Gauss-Hermite moments for velocity profiles derived from STIS 
data. Radii are given in arcsec, first and second moments are given in imits of 
km s~'. 



TABLE 8 
Velocity Profile for NGC 3945 



Radius 


V 






o 


hi 


hi 


0.00 


16.9±10 


.1 


170, 


.7±8.1 


0.032±0.036 


-0.073±0.030 


-0.05 


-14.0±11, 


.4 


173 


.3±9.8 


0.058±0.035 


-0.100±0.029 


-0.10 


-36.0±16, 


.9 


164, 


.4±16.5 


0.146±0.044 


-0.002±0.045 


-0.18 


-46.6±16 


.2 


148, 


.0±18.9 


0.111±0.059 


-0.015±0.057 


-0.30 


-24.0±18, 


.5 


125, 


.4±22.5 


0.093±0.058 


-0.036±0.062 


-0.58 


-100.1±19, 


.0 


138, 


.9±24.2 


-0.038±0.073 


-0.037±0.079 


-1.12 


-29.4±20, 


.5 


138, 


.5±20.5 


-0.035±0.057 


-0.063±0.075 


0.05 


39.5±10, 


.4 


171 


.1±9.3 


0.000±0.034 


-0.074±0.024 


0.10 


12.8±11, 


.8 


158, 


.6±9.9 


-0.021±0.039 


-0.075±0.036 


0.18 


49.0±13 


.1 


121, 


.3±10.9 


-0.011±0.043 


-0.077±0.033 


0.30 


61.3±18, 


.0 


140 


.5±18.0 


-0.053±0.069 


-0.025±0.048 


0.58 


67.1±19 


,7 


167 


.0±17.7 


-0.052±0.059 


-0.077±0.064 


1.12 


101.0±16 


.1 


153, 


.6±16.1 


0.062±0.064 


-0.009±0.055 



Note. — Gauss-Hermite moments for velocity profiles derived from STIS 
data. Radii are given in arcsec, first and second moments are given in units of 
km s~^ 
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TABLE 9 
Velocity Prohle for NGC 4026 



Radius 


V 






a 


hi 


h4 


0.00 


-52.0±23, 


.1 


258 


.3±32.0 


-0.019±0.066 


0.139±0.078 


0.05 


-60.0±26, 


.7 


197, 


.9±27.0 


-0.033±0.059 


-0.001±0.050 


0.10 


-75.6±27, 


.3 


197, 


.1±24.2 


0.098±0.050 


-0.004±0.050 


0.18 


-100.6±27, 


.0 


178, 


.8±30.6 


0.069±0.059 


0.010±0.060 


0.30 


-79.7±27, 


.5 


163, 


.3±22.0 


0.048±0.049 


-0.057±0.053 


0.58 


-123.2±23, 


.2 


176, 


.8±21.6 


0.018±0.056 


-0.066±0.029 


1.12 


-105.5±29, 


.8 


187, 


.5±31.6 


0.()00±0.061 


-0.074±0.060 


-0.05 


-31.5±36, 


.5 


360, 


.5±42.1 


0.158±0.102 


0.099±0.104 


-0.10 


34.5±27, 


.1 


237, 


.7±28.4 


-0.004±0.066 


0.010±0.063 


-0.18 


148.2±24, 


.1 


173, 


.7±26.3 


0.002±0.059 


-0.049±0.051 


-0.30 


189.8±20, 


.7 


155, 


.6±16.4 


-0.006±0.034 


-0.052±0.023 


-0.58 


170.9±18 


.5 


143, 


.9±16.7 


0.007±0.049 


-0.044±0.020 


-1.12 


186.4±15, 


.3 


132 


.1±14.4 


-0.041±0.037 


-0.040±0.034 



Note. — Gauss-Hermite moments for velocity profiles derived from STIS 
data. Radii are given in arcsec, first and second moments are given in units of 

km s"'. 



TABLE 10 





Vhloc 


'1 TY PROl-ILli 


l-OR NGC 5576 




Radius 


V 


O 


h3 




0.00 
-0.05 
-0.10 
-0.18 
-0.30 

0.05 


-1.2±17.3 
17.7±15.7 
-66.7±25.4 
-38.0±18.5 
-17.6±20.0 

-1.5±18.7 


226.4±18.3 
217.7±22.6 
183.6±85.9 
202.1 ±22.6 
218.8±15.1 
215.8±20.9 


0.088±0.054 
0.057±0.051 
-0.016±0.147 
-0.117±0.049 
-0.156±0.048 
0.000±0.000 


0.059±0.054 
0.039±0.048 
-0.035±0.173 
0.041 ±0.034 
0.019±0.042 
0.000±0.000 



Note. — Gauss-Hermite moments for velocity profiles derived from STIS 
data. Radii are given in arcsec, first and second moments are given in imits of 
km s~'. 



TABLE 1 1 

Ground-Based Velocity Profile for NGC 3945 



Axis 


Radius 




V 




a 


hi 


h4 




Major 


0.00 


0, 


,0±0.0 


199 


.9±2.7 


0.016±0.003 


0.147±0. 


,147 


Major 


0.37 


19, 


,3±4.6 


176 


.4±2.5 


0.003±0.010 


0.068±0. 


,109 


Major 


0.74 


42, 


,5±11.0 


180 


.9±3.5 


-0.051±0.010 


0.065±0. 


,096 


Major 


1.30 


52, 


,9±13.3 


180 


.5±5.8 


-0.052±0.026 


0.070±0. 


,083 


Major 


2.04 


75, 


,8±12.5 


201 


.2±6.9 


-0.054±0.049 


0.093±0, 


115 


Major 


3.15 


119, 


,2±14.0 


176, 


.0±14.5 


-0.062±0.044 


0.184±a 


199 


Major 


5.01 


154, 


,1±9.2 


175 


.9±7.9 


-0.121±0.014 


0.215±a 


,223 


Major 


7.98 


163, 


,3±21.5 


286 


.5±122.9 


-0.071±0.036 


o.3ii±a 


,330 


Minor 


0.00 


0, 


,0±0.0 


183 


.9±1.7 


-0.059±0.006 


0.043±0, 


,043 


Minor 


0.37 


21, 


,7±14.1 


176 


.3±8.4 


-0.042±0.036 


0.070±0, 


,092 


Minor 


0.74 


-2, 


,7±1.7 


193 


.5±6.0 


-0.025±0.036 


0.093±0, 


106 


Minor 


1.30 


2, 


,1±4.6 


178 


.3±5.4 


-().()27±0.049 


0.073±0, 


,092 


Minor 


2.04 


-2, 


,4±17.8 


161, 


.7±2.2 


0.004±0.072 


0.089±0, 


,090 


Minor 


3.15 


-6, 


,3±11.4 


142, 


.9±6.4 


-0.024±0.059 


O.I80±0, 


199 



Note. — Gauss-Hermite moments for velocity profiles derived from MDM data. 
Radii are given in arcsec, first and second moments are given in units of km s~^. 
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TABLE 12 

Ground-Based Velocity Profile for NGC 5576 



Axis 


Radius 


y 




a 


hi 


/i4 


Major 


0.00 


3.1±3.5 


188, 


,8±4.6 


0.027±0.021 


-0.023±0.028 


Major 


0.25 


-1.0±3.2 


192, 


,9±4.9 


-0.018±0.020 


-0.015±0.020 


Major 


0.50 


2.1±3.3 


198, 


,4±5.3 


-0.012±0.022 


-0.007±0.026 


Major 


0.88 


4.3±3.7 


191 


,3±5.8 


0.013±0.020 


-0.004±0.023 


Major 


1.38 


8.9±4.6 


184, 


,6±5.6 


-0.080±0.021 


-0.013±0.027 


Major 


2.12 


10.0±4.5 


166, 


,8±7.9 


-0.045±0.023 


0.118±0.022 


Major 


3.38 


28.9±5.4 


166, 


,3±8.9 


-0.113±0.030 


0.060±0.035 


Major 


5.38 


16.7±7.8 


162, 


,6±12.3 


-0.041±0.053 


0.006±0.052 


60° Skew 


0.00 


2.0±3.6 


191 


,6±5.6 


-0.029±0.021 


0.011±0.025 


60° Skew 


0.25 


4.9±3.4 


188, 


,9±5.5 


-0.017±0.023 


-0.009±0.023 


60° Skew 


0.50 


3.3±3.8 


188, 


.2±5.1 


-0.012±0.019 


-0.014±0.026 


60° Skew 


0.88 


2.6±3.7 


191, 


,2±5.9 


-0.016±0.023 


0.010±0.027 


60° Skew 


1.38 


3.3±4.5 


192, 


,0±6.7 


-0.012±0.023 


0.018±0.025 


60° Skew 


2.12 


9.6±4.1 


190, 


.1±4.5 


-0.008±0.026 


-0.062±0.024 


60° Skew 


3.38 


11.0±5.0 


171, 


.7±5.6 


-0.102±0.027 


-0.010±0.035 


60° Skew 


5.38 


23.1±5.3 


146 


,7±7.0 


-0.110±0.028 


0.021 ±0.032 


60° Skew 


-0.25 


3.1±3.8 


195, 


,1±6.3 


0.026±0.020 


0.025±0.026 


60° Skew 


-0.50 


1.0±4.2 


196, 


,0±8.4 


0.036±0.021 


0.050±0.027 


60° Skew 


-0.88 


1.4±4.6 


183, 


,8±6.7 


0.051 ±0.022 


0.027±{).{)25 


60° Skew 


-1.38 


15.7±4.0 


168, 


,1±6.2 


0.032±0.023 


-0.018±0.023 


60° Skew 


-2.12 


7.4±4.4 


159, 


,6±6.6 


0.025±0.028 


-0.026±0.030 


60° Skew 


-3.38 


10.5±5.0 


158, 


,2±7.9 


0.048±0.028 


0.013±0.028 


60° Skew 


-5.38 


13.8±6.2 


157, 


.5±7.7 


-0.028±0.034 


0.018±0.029 



Note. — Gauss-Hermite moments for velocity profiles derived from Magellan 
data. Radii are given in arcsec, first and second moments are given in units of km s~^ 



TABLE 13 

Ground-Based Surface Brightness Prohle 
for NGC 3945 



Radius 


Surface Brightness 


Ellipticiy 


RA. 


10.28 


18.909 


0.352 


156.1 


12.10 


19.226 


0.328 


155.8 


14.23 


19.664 


0.236 


154.8 


16.75 


20.052 


0.140 


154.7 


19.70 


20.416 


0.038 


94.4 


23.18 


20.591 


0.134 


73.7 


27.27 


20.678 


0.262 


70.6 


32.08 


20.929 


0.300 


70.1 


37.74 


21.517 


0.175 


74.6 


44.40 


21.995 


0.100 


145.3 


52.23 


22.415 


0.169 


150.0 


61.45 


22.928 


0.089 


164.3 



Note. — Radius is given in units of arcsec. Sur- 
face brightness is V-band in units of magnitudes per 
square arcsec. The third column gives ellipticity, and 
the fourth column gives position angle in degrees east 
of north. 



